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I. HEAVY QUARK SYMMETRY 



In these lectures I will introduce the ideas of heavy quark symmetry and the heavy quark 
limit, which exploit the simplification of certain aspects of QCD for infinite quark mass, 
rriQ — > oo. We will see that while these ideas are extraordinarily simple from a physical 
point of view, they are of enormous practical utility in the study of the phenomenology of 
bottom and charmed hadrons. One reason for this is the existence not just of an interesting 
new limit of QCD, but of a systematic expansion about this limit. The technology of this 
expansion is the Heavy Quark Effective Theory (HQET), which allows one to use heavy quark 
symmetry to make accurate predictions of the properties and behavior of heavy hadrons in 
which the theoretical errors are under control. While the emphasis in these lectures will be 
on the physical picture of heavy hadrons which emerges in the heavy quark limit, it will 
be important to introduce enough of the formalism of the HQET to reveal the structure 
of the heavy quark expansion as a simultaneous expansion in powers of Aqcd/^q and 
a s (mg). However, what I hope to leave you with above all is an appreciation for the 
simplicity, elegance and coherence of the ideas which underlie the technical results which 
will be presented. The interested reader is also encouraged to consult a number of excellent 
reviews |T| , which typically cover in more detail the material in the first two sections of these 
lectures. 



Why is an understanding of QCD crucial to the study of the properties of the bottom 
and charmed quarks? As an example, consider semileptonic b decay, b — ► civ. This process 
is mediated by a four-fermion operator, 



A. Introduction 




(1.1) 



The weak matrix element is easy to calculate at the quark level, 
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Ajuark = (du\ O bc \b) = —^U c (p c )^(l - rf)u b {p b ) U^p^^l - ^)v u (p v ) . (1.2) 

However, ^4 qua rk is only relevant at very short distances; at longer distances, QCD confine- 
ment implies that free b and c quarks are not asymptotic states of the theory. Instead, 
nonperturbative QCD effects "dress" the quark level transition b — > civ to a, hadronic 
transition, such as 

B -> Dlv or B -> D*tv or ... (1.3) 

(In these lectures, we will use a convention in which a B meson contains a b quark, not a 
b antiquark.) The hadronic matrix element w4h a dron depends on nonperturbative QCD as 
well as on GpV c b, and is difficult to calculate from first principles. To disentangle the weak 
interaction part of this complicated process requires us to develop some understanding of 
the strong interaction effects. 

There are a variety of methods by which one can do this. Perhaps the most popular, 
historically, has been use of various quark potential models @. While these models are 
typically very predictive, they are based on uncontrolled assumptions and approximations, 
and it is virtually impossible to estimate the theoretical errors associated with their use. 
This is a serious defect if one builds such a model into the experimental extraction of a weak 
coupling constant such as V c b, because the uncontrolled theoretical errors then infect the 
experimental result. 

Another approach, the one to be discussed here, is to exploit as much as possible the fact 
that the b and c quarks are heavy, by which we mean that m&, m c 3> Aqcd- The scale Aqcd 
is the typical energy at which QCD becomes nonperturbative, and is of the order of hundreds 
of MeV. The physical quark masses are approximately m& pa 4.8 GeV and m c pa 1.5 GeV. 
The formalism which we will develop will not make as many predictions as do potential 
models. However, the compensation will be that we will develop a systematic expansion 
in powers of Aqcd/^6,c, within which we will be able to do concrete error analysis. In 
particular, we will be able to estimate the error associated with the fact that m c may not be 
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very close to the asymptotic limit m c 3> Aqcd- Even where this error may be substantial, 
the fact that it is under control allows us to maintain predictive power in the theory. 

B. The heavy quark limit 

Consider a hadron Hq composed of a heavy quark Q and "light degrees of freedom", 
consisting of light quarks, light antiquarks and gluons, in the limit mg — > oo. The Compton 
wavelength of the heavy quark scales as the inverse of the heavy quark mass, Aq ~ 1/mq. 
The light degrees of freedom, by contrast, are characterized by momenta of order Aqcd; 
corresponding to wavelengths A^ ~ 1/Aq CD . Since A^ Aq, the light degrees of freedom 
cannot resolve features of the heavy quark other than its conserved gauge quantum numbers. 
In particular, they cannot probe the actual value of Aq, that is, the value of mq. 

We draw the same conclusion in momentum space. The structure of the hadron Hq is 
determined by nonperturbative strong interactions. The asymptotic freedom of QCD implies 
that when quarks and gluons exchange momenta p much larger than Aqcd> the process is 
perturbative in the strong coupling constant a s (p). On the other hand, the typical momenta 
exchanged by the light degrees of freedom with each other and with the heavy quark are of 
order Aqcd, for which a perturbative expansion is of no use. For these exchanges, however, 
p < mq, and the heavy quark Q does not recoil, remaining at rest in the rest frame of the 
hadron. In this limit, Q acts as a static source of electric and chromoelectric gauge field. 
The chromoelectric field, which holds Hq together, is nonperturbative in nature, but it is 
independent of mq. The result is that the properties of the light degrees of freedom depend 
only on the presence of the static gauge field, independent of the flavor and mass of the 
heavy quark carrying the gauge charge.[] 

There is an immediate implication for the spectroscopy of heavy hadrons. Since the 



Top quarks decay too quickly for a static chromoelectric field to be established around them, so 
the simplifications discussed here are not relevant to them. 
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FIG. 1. Schematic spectra of the bottom and charmed mesons in the limit mb,m c 3> Aqcd- 
The offset of the two spectra is not to scale; in reality, m& — m c 3> Aj ~ Aqcd- 

interaction of the light degrees of freedom with the heavy quark is independent of tuq, then 
so is the spectrum of their excitations. It is these excitations which determine the spectrum 
of heavy hadrons Hq. Hence the splittings Ai ~ Aqcd between the various hadrons Hq are 
independent of Q and, in the limit rriQ — > oo, do not scale with tuq. 

For example, the spectra of bottom and charmed mesons are shown schematically in 
Fig. [I], in the limit mf,,m c 3> Aqcd- The light degrees of freedom are in exactly the same 
state in the mesons B{ and Di, for a given i. The offset Bi — Di = nib — ™>c is just the 
difference between the heavy quark masses; in no way does the relationship between the 
spectra rely on an approximation m& ~ m c . 

We can enrich this picture by recalling that the heavy quarks and light degrees of freedom 
also carry spin. The heavy quark has spin quantum number Sq = ~, which leads to a 
chromomagnetic moment 



HQ 



2rrif 



Note that [Xq —> as mg — > oo, and the interaction between the spin of the heavy quark and 
the light degrees of freedom is suppressed. Hence the light degrees of freedom are insensitive 



to Sq; their state is independent of whether Sq = \ or Sq = — ~. Thus each of the energy 
levels in Fig. [I] is actually doubled, one state for each possible value of Sq. 

To summarize, what we see is that the light degrees of freedom are the same when 
combined with any of the following heavy quark states: 

Qi(T), 2 (t), a (i); ••• Q* h (T), Qw k U), (1.5) 

where there are iV^ heavy quarks (in the real world, Nh = 2). The result is an SU(2Nh) 
symmetry which applies to the light degrees of freedom 0-0. A new symmetry means new 
nonperturbative relations between physical quantities. It is these relations which we wish 
to understand and exploit. 

The light degrees of freedom have total angular momentum Jt, which is integral for 
baryons and half-integral for mesons. When combined with the heavy quark spin Sq = |, 
we find physical hadron states with total angular momentum 



J 



;i.6) 



If St 7^ 0, then these two states are degenerate. For example, the lightest heavy mesons 
have Si = |, leading to a doublet with J = and J — 1. Indeed, in the charm system we 
find that the states of lowest mass are the spin-0 D and the spin-1 D*; the corresponding 
bottom mesons are the B and B*. 

When effects of order l/rriQ are included, the chromomagnetic interactions split the 
states of given St but different J. This "hyperfine" splitting is not calculable perturbatively, 
but it is proportional to the heavy quark magnetic moment /xq. This gives its scaling with 

mo* — vriD ~ l/m c 

m B * -m B ~ l/m b . (1.7) 

From this fact we can construct a relation which is a nonperturbative prediction of heavy 
quark symmetry, 



m 2 B * — m 2 B = m 2 D « — m 2 D . (1.8) 

Experimentally, m B * — m 2 B = 0.49 GeV 2 and m 2 D * — m 2 D = 0.55 GeV 2 , so this prediction 
works quite well. Note that this relation involves not just the heavy quark symmetry, but 
the systematic inclusion of the leading symmetry violating effects. 

So far, we have formulated heavy quark symmetry for hadrons in their rest frame. Of 
course, we can easily boost to a frame in which the hadrons have arbitrary four-velocity 
yt 1 = ry(\ 7 For heavy quarks Q\ and Q2, the symmetry will then relate hadrons H\(v) and 
H 2 {v) with the same velocity but with different momenta. This distinguishes heavy quark 
symmetry from ordinary symmetries of QCD, which relate states of the same momentum. 
To remind ourselves of this distinction, henceforth we will label heavy hadrons explicitly by 
their velocity: D(v), D*(v), B(v), B*(v), and so on. 

C. Semileptonic decay of a heavy quark 

Now let us return to the semileptonic weak decay b — > civ, but now consider it in the 
heavy quark limit for the b and c quarks. Suppose the decay occurs at time t = 0. For t < 0, 
the b quark is embedded in a hadron Hj,; for t > 0, the c quark is dressed by light degrees 
of freedom to H' c . Let us consider the lightest hadrons, H\, = B{y) and H' c = D{v'). Note 
that since the leptons carry away energy and momentum, in general v 7^ v'. 

What happens to the light degrees of freedom when the heavy quark decays? For t < 0, 
they see the chromoelectric field of a point source with velocity v. At t — 0, this point 
source recoils instantaneously^ to velocity v '; the color neutral leptons do not interact with 
the light hadronic degrees of freedom as they fly off. The light quarks and gluons then 
must reassemble themselves about the recoiling color source. This nonperturbative process 
will generally involve the production of an excited state or of additional particles; the light 
degrees of freedom can exchange energy with the heavy quark, so there is no kinematic 



2 The weak decay occurs over a very short time 5t ~ 1/Mw <C 1/Aqcd- 
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restriction on the excitations (of energy ~ Aqcd) which can be formed. There is also some 
chance that the light degrees of freedom will reassemble themselves back into a ground state 
D meson. The amplitude for this to happen is a function only of the inner product w = v ■ v ' 
of the initial and final velocities of the color sources. This amplitude, C,(w), is known as the 
Isgur-Wise function [^J. 

Clearly, the kinematic point v = v', or w = 1, is a special one. In this corner of phase 
space, where the leptons are emitted back to back, there is no recoil of the source of color 
field at t = 0. As far as the light degrees of freedom are concerned, nothing happens! Their 
state is unaffected by the decay of the heavy quark; they don't even notice it. Hence the 
amplitude for them to remain in the ground state is exactly unity. This is reflected in a 
nonperturbative normalization of the Isgur-Wise function at zero recoil [[3], 

£(1) = 1. (1.9) 

As we will see, this normalization condition is of enormous phenomenological use. It will 
be extremely important to understand the corrections to this result for finite heavy quark 
masses mj and, especially, m c . 

The weak decay b — > c is mediated by a left-handed current C7 M (1 — / y 5 )b. Not only 
does this operator carry momentum, but it can change the orientation of the spin Sq of 
the heavy quark during the decay. For a fixed light angular momentum Jg, the relative 
orientation of Sq determines whether the physical hadron in the final state is a D or a D* . 
However, the light degrees of freedom are insensitive to Sq, so the nonperturbative part 
of the transition is the same whether it is a D or a D* which is produced. Hence heavy 
quark symmetry implies relations between the hadronic matrix elements which describe the 
semileptonic decays B — > Div and B — ► D*lv. 

It is conventional to parameterize these matrix elements by a set of scalar form factors. 
These are defined separately for the vector and axial currents, as follows: 

(D(v')\ cry»b \B(v)) = h+(w)(v + vj l + h-(w)(v - v'Y 
(D*(v',e)\ cYb\B(v)) = h v (w)ie^elv' a v p 
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(D(v')\cy»7 5 b\B(v)) = (1.10) 
(D(v', e)\cYl 5 b \B(v)) = h Al {w){w + l)e*" - e* ■ v[h Aa (w)v» + h A3 {w)v'»} . 

The set of form factors hi(w) is the one appropriate to the heavy quark limit. Other linear 
combinations are also found in the literature. In any case, the form factors are independent 
nonperturbative functions of the recoil or equivalently, for fixed mj and m c , of the momentum 
transfer. However, in the heavy quark limit they correspond to a single transition of the 
light degrees of freedom, being distinguished from each other only by the relative orientation 
of the spin of the heavy quark. Hence they may all be written in terms of the single function 
which describes this nonperturbative transition. As we will derive later, the result is a 
set of relations ||, 

h+(w) = h v (w) = h Al (w) = h A3 (w) = £(w) 

h_(w) = h A2 (w) = 0, (1.11) 

which follow solely from the heavy quark symmetry. Of course, all of the form factors 
which do not vanish inherit the normalization condition ( |1.9| ) at zero recoil. This result is a 
powerful constraint on the structure of semileptonic decay in the heavy quark limit. 



D. Heavy meson decay constant 

As a final example of the utility of the heavy quark limit, consider the coupling of the 
heavy meson field to the axial vector current. This is conventionally parameterized in terms 
of a decay constant; for example, for the B~ meson we define f B via 

(0\uYl 5 b\B-(p B )) = if B p%. (1.12) 

What is the dependence of the nonperturbative quantity f B on m B 7 To address this question, 
we rewrite Eq. ( |1 . 1 2| ) in a form appropriate to taking the heavy quark limit, m B — > oo (which 
is equivalent to m;, — > oo). This entails making explicit the dependence of all quantities on 
m B . First, we trade the B~ momentum for its velocity, 
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p B = m B v^. (1.13) 
Second, we replace the usual B~ state, whose normalization depends on m B , 

(B(p 1 )\B(p 2 )) = 2E B 5 (3 \p 1 -p 2 ), (1.14) 
by a mass-independent state, 

\B{v)) = -±=\B(p B )), (1.15) 

satisfying 

(B( Vl )|S(v 2 )) = 2 7 5( 3 )(p 1 -p 2 ). (1.16) 

Then Eq. ( |1.12| ) becomes 

v^(0| u 7 V& = if B m B v» . (1.17) 

The nonperturbative matrix element (0| icy^^b \B~(v )) is independent of mg in the heavy 
quark limit. Hence, we see that in this limit f B takes the form 

f B = m B x (independent of m B ) . (1-18) 

This makes explicit the scaling of f B with m B . It is more interesting to write this as a 
prediction for the ratio of charmed and bottom meson decay constants. We find 



fB _ /^ +0 MqcD ; AqCd\ (119) 



f D y m B V m D m B 
For the physical bottom and charm masses, of course, the correction terms proportional to 
Aqcd/^iq could be important. 

II. HEAVY QUARK EFFECTIVE THEORY 

We have already extracted quite a bit of nontrivial information from the heavy quark 
limit. We have found the scaling of various quantities with m,Q, we have studied the impli- 
cations for heavy hadron spectroscopy, and we have found nonperturbative relations among 
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the hadronic form factors which describe semileptonic b — > c decay. However, all of these 
results have been obtained in the strict limit mq —>■ oo. If the heavy quark limit is to be 
of more than academic interest, and is to provide the basis for quantitative phenomenol- 
ogy, we have to understand how to include corrections systematically. There are actually 
two types of corrections which we would like to include. Power corrections are subleading 
terms in the expansion in Aqcd/^q; those proportional to Aqcd/^c are the most wor- 
risome, because of the relatively small charm quark mass. Logarithmic corrections arise 
from the implicit dependence of quantities on mq through the strong coupling constant 
« s (wq) ~ 1/ ln(mg/AQCD)- For the physical values of m;, and m c , either of these could be 
important. What we need is a formalism which can accommodate them both. 

In short, we need to go from a set of heavy quark symmetry predictions in the mq — > oo 
limit, to a reformulation of QCD which provides a controlled expansion about this limit. 
The formalism which does the job is the Heavy Quark Effective Theory, or the HQET. The 
purpose of the HQET is to allow us to extract, explicitly and systematically, all dependence 
of physical quantities on mq, in the limit mq ^> Aqcd- In these lectures, we will develop 
only enough of the technology to treat the dominant leading effects, providing indications 
along the way of how one would carry the expansion further. 

The HQET, as formulated here, was developed in a series of papers going back to the late 
1980's |3|.|^J^,3-[T5| , which the reader who is interested in tracking its historical development 
may consult. 



Consider the kinematics of a heavy quark Q, bound in a hadron with light degrees of 
freedom to make a color singlet state. The small momenta which Q typically exchanges with 
the rest of the hadron are of order Aqcd <C m Qi and they never take Q far from its mass 



A. The effective Lagrangian 



shell, pq = m, 



2 q. Hence the momentum pq can be decomposed into two parts, 



Pq = m Q vl1 + ^ 



(2.1) 
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where rnQV^ is the constant on-shell part of pg, and k^ ~ Aqcd is the small, fluctuating 
"residual momentum" . The on-shell condition for the heavy quark then becomes 

m 2 Q = (m Q v fl + k^f = m 2 Q + 2m Q v ■ k + k 2 . (2.2) 

In the heavy quark limit we may neglect the last term, and we have the simple condition 

v ■ k = (2.3) 

for an on-shell heavy quark. Here the velocity functions as a label; since soft interactions 
cannot change v^, there is a velocity superselection rule in the heavy quark limit, and is 
a good quantum number of the QCD Hamiltonian. 

We find the same result by taking the mg — > oo limit of the heavy quark propagator, 

i 1 + j) i 



p 1 — rriQ + ie 2 v ■ k + ie 
In this limit the propagator is independent of toq. The projection operators 



(2.4) 



project onto the positive (P + ) and negative (P_) frequency parts of the Dirac field Q. 
This is clear in the Dirac representation in the rest frame, in which P + and P_ project, 
respectively, onto the upper two and lower two components of the heavy quark spinor. In 
the limit uiq — > oo, in which Q remains almost on shell, only the "large" upper components 
of the field Q propagate; mixing via zitterbewegung with the "small" lower components is 
suppressed by l/2rriQ. Hence the action of the projectors on Q is 

P + Q(x) = Q(x) + 0(l/m Q ), P-Q(x) = + 0(l/m Q ) . (2.6) 

The momentum dependence of the field Q is given by its action on a heavy quark state, 

Q(x) \Q(p)} = e-** |0> . (2.7) 

If we now multiply both sides by a phase corresponding to the on-shell momentum, 
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e im ^ x Q(x) \Q(p)) = e~ ik - x |0) , (2.8) 

the right side of this equation is independent of mq. Hence the left side must be, as well. 
Combining this observation with the argument of the previous paragraph, we are motivated 
to define a mq- independent effective heavy quark field h v (x), 

h v (x) = e tm ^ x P + Q(x). (2.9) 

Note that the effective field carries a velocity label v and is a two-component object. The 
modifications to the ordinary field Q(x) project out the positive frequency part and ensure 
that states annihilated by h v (x) have no dependence on mq. Hence, these are reasonable 
candidate fields to carry representations of the heavy quark symmetry. Of course, the small 
components cannot be neglected when effects of order 1/mq are included. In the HQET 
they are represented by a field 

H v (x) = e im ^ x P^Q{x). (2.10) 

The field H v (x) vanishes in the mq — > oo limit. 

The ordinary QCD Lagrange density for a field Q(x) is given by 

£qcd = Q(x) (ip - m Q ) Q(x) , (2.11) 

where = — igA®T a is the gauge covariant derivative. To find the Lagrangian of the 
HQET, we substitute 

Q(x) = e- im ^ x h v {x) + . . . (2.12) 
into £qcd and expand. With the aid of the projection identity P + ^P + = v^, we find 

£ H qet = h v (x)iv ■ Dh(x) . (2.13) 
This simple Lagrangian leads to the propagator we derived earlier, 

tot < 2 ' 14 > 
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and to an equally simple quark-gluon vertex, 

igTVA^. (2.15) 

Note that both the propagator and the vertex are independent of ijiq, reflecting the heavy 
quark flavor symmetry. They also have no Dirac structure, reflecting the heavy quark 
spin symmetry. Our intuitive statements about the structure of heavy hadrons have been 
promoted to explicit symmetries of the QCD Lagrangian in the limit mg — > oo. 

It is straightforward to include power corrections to £hqet- Write Q(x) in terms of the 
effective fields, 

Q( x ) = e~ m ^ x [h v (x) + H v (x)} , (2.16) 
and apply the classical equation of motion (ip — mQ)Q[x) = 0: 

iph v {x) + (ip -2m Q )H v (x) = 0. (2.17) 
Multiplying by P_ and commuting f to the right, we find 

(iv ■ D + 2m Q )H v (x) = ip± h v (x) , (2.18) 

where = — v^v ■ D. We then substitute Q(x) into £qcd as before, eliminate H v (x) 
and expand in 1/mq to obtain 

^hqet = Kiv ■ Dh v + h v ip± — ^— ip±h v 

iv ■ D + 2rriQ 



h v iv ■ Dh v + 



2mQ 



h v (iD ± ) 2 h v + ^-h v a al3 G a 8h v 



+ .... (2.19) 



2 

The leading corrections have a simple interpretation, which becomes clear in the rest frame, 
f M = (1,0,0,0). The spin-independent term is 

7^—O k = -^— h v (iD ± ) 2 h v -> h v {iDfh v , (2.20) 

2m Q 2m Q 2m Q 

which is the negative of the nonrelativistic kinetic energy of the heavy quark. Because of the 
explicit factor of l/2m.Q, this term violates the heavy flavor symmetry. The spin-dependent 
part is 
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O g = —\ h v a^G aP h v - h v a ij T a h v x gG% = gfl a Q ■ B a , (2.21) 

2m Q 2m Q 2 Am Q J v 

which is the coupling of the spin of the heavy quark to the chromomagnetic field. Because 
it has a nontrivial Dirac structure, this term violates both the heavy flavor symmetry and 
the heavy spin symmetry. For example, Oq is responsible for the D — D* and B — B* mass 
splittings. These correction terms will be treated as part of the interaction Lagrangian, even 
though Ok has a piece which is a pure bilinear in the heavy quark field. 



B. Effective currents and states 

The expansion of the weak interaction current C7 M (1 — 7 5 )6 is analogous. However, here 
we must introduce separate effective fields for the charm and bottom quarks, each with its 
own velocity: 

b^h b v , c^h c v ,. (2.22) 

Then a general flavor- changing current becomes, to leading order, 

cTb^h c v ,Th b v} (2.23) 

where T is a fixed Dirac structure. With the leading power corrections, this is 

cT b - hi, Th b v + J- K,T(ip ± )h b v + ^- \)Yh b v + .... (2.24) 

The effective currents, and other operators which appear in the HQET, may often be sim- 
plified by use of the classical equation of motion, 

iv-Dh v (x)=0. (2.25) 

However, it is only safe to apply these equations naively at order 1/mg; at higher order the 



application of the equations of motion involves additional subtleties p)|-|T8|j . 

To complete the effective theory, we need mg-independent hadron states which are cre- 
ated and annihilated by currents containing the effective fields. For example, there is an ef- 
fective pseudoscalar meson state \M(v)) which couples to the effective axial current qj^j^hv, 
with a coupling Fm which is independent of tuq, 
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FIG. 2. Tree level plus the one loop renormalization of the current qTQ in QCD. The box 

represents the current insertion. 



(0\q~f' t 'fh v \M(v)) = iF M v». 



(2.26) 



At lowest order, Fm is related to a conventional decay constant such as fs by 



Sb = F M /\/rnB~ 



(2.27) 



from which we immediately find the relationship (|1.19|) between fjj and 



C. Radiative corrections 

We can use the effective Lagrangian £hqet to compute the radiative corrections to 
the matrix element ( j2.26j ). In particular, we would like to extract the dependence of F M on 
In m,Q. This dependence comes through the one-loop renormalization of the current q^^^Q. 
At lowest order, of course, the renormalization is straightforward: we simply compute the set 
of graphs found in Fig. |[ The result is finite, because the current is (partially) conserved, 
and we extract a result of the form 

qri 5 Q (l - 7o^ Hm Q /m q ) + ...). (2.28) 

Note that there is no explicit dependence on the renormalization scale fi, since there is no 
divergence to be subtracted. 

The same result may be obtained in the effective theory. In this case we must match 
the currents in full QCD onto HQET currents of the form q r )^ r j 5 h v . This step will induce a 
matching coefficient containing the explicit dependence on ttlq, which is absent, by construc- 
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h, q h, q h, q h, <l 

FIG. 3. Tree level plus the one loop renormalization of the current qj^^hy in HQET. The 

double line is the propagator of the effective field h v . 

tion, from the operators and Lagrangian of the HQET.Q In addition, the effective current 
will not necessarily be conserved, since the ultraviolet properties of QCD and the HQET 
differ. Hence the form of the matching, once radiative corrections are included, is 

qn 5 Q - C(m Q , /i) qn'K . (2.29) 

We can deduce the form of C(rriQ,fjj) by considering the renormalization of the effective 
current q'y^^hy, shown in the last three terms of Fig. |3|. These diagrams, computed in the 
effective theory, are independent of rriQ. However, in general they are divergent, so they 
depend on the renormalization scale \x\ the renormalization takes the form 

qYl 5 K (c(m Q , fi) - 7o ^ lnfa/m,) + ...), (2.30) 

where here m q acts as an infrared cutoff. The \x dependence in the second term must be 
canceled by C(uiq, h). Since the logarithm depends on a dimensionless ratio, C(rriQ, fi) must 
be of the form 

a 

C{m Q , fi) = 1- 7o— ln(m Q /fj,) + ... . (2.31) 

7T 

Comparing the dependence on uiq of C(rriQ, fi) and the expansion (|2.28| ), we see that 7 = 7o- 
However, the effective theory allows us to go beyond leading order, and to resum all 
corrections of the form a™ln n mg. We do this with the renormalization group equations, 

3 In general, the matching procedure at order a s can also induce new Dirac structures qT h v . They 
do not affect the leading logarithms discussed here. 
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which express the independence of physical observables on the renormalization scale /i. In 
this case, they require that the /i dependence of C(mg,/i) cancel that of the one loop 
diagrams in Fig. ^ under small changes in /i: 

/^C(m Q ,/,)= 7 o^. (2.32) 

Here the total derivative with respect to // includes the implicit dependence on \i of the 
coupling constant a s {\i): 

d d n , . d 

2 2^ 

/3 = 11 - = — foriV / = 4, (2.33) 

where Nf is the number of light flavors. We compute the anomalous dimension 70 from the 
divergent parts of the one loop diagrams shown in Fig. |3|. From the three terms we find, 
respectively H||, 

8 1 / 8\ 1 /16\ 

70 = 77+0 -o +ohr = 4 - ( 2 - 34 ) 



3 2 V 3/ 2 V 3 
The solution of the renormalization group equation is 

,(m Q )y o/2 "° = ( a s (m Q ) 



This then yields the leading logarithmic correction to the ratio 

6/25 



/b \m D ( a s (m b ) 



(2.36) 



/d V m B \oc s {m c ) J 

The radiative correction is approximately a ten percent effect. In fact, it has a simple 
physical interpretation. For virtual gluons of "intermediate" energy, m c < E g < nib, the 
bottom quark is heavy but the charm quark is light. Such gluons contribute to the difference 
between and fjj even in the heavy quark limit. 

In summary, then, the purpose of the HQET is to make explicit all dependence of ob- 
servable quantities on tuq. The logarithmic dependence, through a s (rriQ), arises from in- 
termediate virtual gluons with m c < E g < m&. We obtain these corrections by computing 



perturbatively with the HQET Lagrangian, then using the renormalization group to resum 
the logarithms to all orders. The power dependence, 1/mg, is extracted systematically in 
the heavy quark expansion. We have seen how to expand the Lagrangian and the states to 
subleading order; the application of the expansion to a physical decay rate will be presented 
in the next section. 

These lectures are meant to be pedagogical, so we will only treat the leading corrections 
to a few processes. However, the state of the art goes significantly beyond what will be 
presented here. For many quantities, not only the leading logarithms, a™ \n n rriQ, but the 
subleading (two loop) logarithms, of order a^ +1 ln n mg, have been resummed. Similarly, 
many power corrections are known to relative order l/rriq. It is particularly important 
phenomenologically to take into account the corrections of order 



III. EXCLUSIVE B DECAYS 

We now have the tools we need for an HQET treatment of the exclusive semileptonic 
transitions B — > D iv and B — > D* tv. Earlier, we argued on physical grounds that in 
the heavy quark limit all of the hadronic matrix elements which appear in these decays 
are related to a single nonperturbative function £(w). Now we will sharpen this analysis 
to actually derive these relations, and to include radiative and power corrections. In fact, 
almost all of our effort will go into the power corrections, since the radiative corrections to 
the transition currents are computed just as in the previous section. 

A. Matrix element relations at leading order 

The transitions in question require the nonperturbative matrix elements 
(D(v')\cYb\B{v)) , (D*(v',e)\&r^b\B(v)) , (£>>', e)| c^^b \B(v)) , (3.1) 



parameterized in terms of form factors as in Eq. ( |1.10| ). Our first task is to derive the relations 



between these form factors, as promised earlier. These relations depend on the heavy quark 
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symmetry, that is, on the fact that the spin quantum numbers of Q and of the light degrees 
of freedom are separately conserved by the soft physics. Hence we need a representation of 
the heavy meson states in which they have well defined transformations separately under 
the angular momentum operators Sq and Je. In particular, the representation must reflect 
the fact that a rotation by Sq can exchange the pseudoscalar meson M(v) with the vector 
meson M*(v, e). 

The solution is to introduce a "superfield" A4(v), defined as the 4x4 Dirac matrix |TTlJl9|| 

1 + 



M(v) 



7^M>,e)-7 & M(t,) . (3.2) 



2 

Under heavy quark spin rotations Sq, M.(v) transforms as 

M(v) -+ D(S Q )M(v) , (3.3) 

and under Lorentz rotations A, as 

M(v) -> D{A)M{A- 1 v)D-\A). (3.4) 

Here D(- • •) is the spinor representation of 5*0(3, 1). The superfield satisfies the matrix 
identity 

P+M(v)P_ = M(v), (3.5) 

so it transforms the same way as the product of spinors h v q, representing a heavy quark 
and a light antiquark moving together at velocity v fl . 

A current which mediates the decay of one heavy quark (Q) into another (Q') is of the 
form h v i T h v . Under a rotation by Sq, the effective field h v transforms as 

h v -> D{S Q ) h v , (3.6) 

while h v i is unchanged. The current would remain invariant if we took V to transform as 

Y^YD-\S Q ). (3.7) 

On the other hand, the matrix element of superfields 

20 



(M'{v')\h V 'Th v \M{v)) (3.8) 

is invariant if we rotate both h v and M.(v) by the same Sq. With the transformation law 
( p.3| ) for A4(v), it follows that the Sg-invariant matrix element must be proportional to 
TA4(v). When we also consider rotations under Sq>, we find that the matrix element is 
restricted to the general form 

(M'(v')\ h v ,Th v \M(v)) = - [M M M M >] 1/2 Tr \m\v')T M(v) F(v,v')] . (3.9) 

The product of masses in front is a convention which restores the relativistic normalization 
of the states. Note that the heavy quark symmetry allows an arbitrary 4x4 Dirac matrix 
F(v, v') to act on the "light quark" part of the superfields. Its presence reflects the fact 
that, other than the constraint of Lorentz symmetry, the behavior of the spin of the light 
degrees of freedom during the decay is unknown. 

A general expansion of F(v, v') in terms of scalar functions Fi(w = v ■ v ') takes the form 

F(v, v') = Fx(w) + F 2 (w)i) + F 3 (w)f + F A (w)ff . (3.10) 

However, the identity (|3.5| ) applied to the matrix element ( ^.9| ) yields 

p(v, v') = P_ F(v, v') PL = [Fi(w) - F 2 {w) - F 3 (w) + F 4 (w)} P_ P'_ . (3.11) 

In other words, F(v, v') is actually a scalar, which we identify with the Isgur-Wise function, 

(3-12) 

As an exercise, let us apply this formalism to the matrix elements for B — > (D, D*) Iv. For a 
given matrix element, we pick out the part of the superfield Ai(v) which is relevant. Hence 
we find 

(D(v')\c^b\B(v)) = (M D (v')\h c v ,Yh b v \M B (v)) 



[m D m B ] 1/2 Tr ^P' + ^P + {-^) i{ 



w 



[m D m B ] l/2 i(w) {v + v'Y (3.13) 
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{D*{v\t)\cn b b\B{v)) = (M* D (v',e)\h c vf n S h h v \M B (v)) 



[m D *m B ] 1/2 Tr 7 5 P + (- 7 5 ) f( 



w 



[m D *m B } 1/2 £(w) [(w + l)e*" - e* • (v + 



(3.14) 



(D*(v', e)| ey M fo |-B(t>)) 



- [m^m B ] 1/2 Tr [/P; 7 ^P+(- 7 5 )] £(w) 



(3.15) 



reproducing explicitly the relations (|1.11| ) between the independent form factors hi{w). We 
can also derive the normalization condition at w = 1. Consider the matrix element of the b 
number current by^b between B meson states. In QCD, the matrix element of this current 
is exactly normalized, 



[B{v)\ by% \B(v)) = 2p B = 2m B v^ 



But in HQET, we have 



(B(v)\brfb \B{v)) = (M B (v)\ h v YK \M b {v)) 
= 2m B v li i(l). 



Hence the normalization condition at zero recoil, 



follows directly from the conservation of the heavy quark number current. 



(3.16) 



(3.17) 



(3.18) 



B. Power corrections to the matrix elements 



The matrix elements we have derived are computed in the strict limit m^^ — > oo. How 
are they affected by corrections of order l/m& and l/m c ? There are two sources of 1/tuq 
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corrections in the effective theory: the corrections ( |2.24|) to the heavy quark currents, and 
the corrections ( |2.19| ) to the Lagrangian. 



When radiative corrections are included, the expansion of the heavy quark current cTb 
in terms of HQET operators has a form which is somewhat more general than Eq. ( |2.24| ), 

cTb - ao(a s ) K,Th b v + K,T1iD a hl + h c v ,{-iD a )T' x a h b v + .... (3.19) 

The matrix elements of the power corrections are constrained by heavy quark symmetry in a 
manner completely analogous to the leading current. In terms of traces over the superfields, 
we have [ |13[ 



(M\v')\h v ,T a iD a h v \M{v)) = -{M M M M ,] 1/2 Tr [m'(v') T a M (v) G a (v, v')] , (3.20) 

where G a (v,v') is another arbitrary 4x4 Dirac matrix. The matrix element 

(M'(v')\ K (-i*D a )T[ a h v \M(v)) (3.21) 

may also be written in terms of G a (v,v'), using charge conjugation. 

The 1/rriQ corrections Ok and Og to the Lagrangian contribute somewhat differently. 
In order to apply heavy quark symmetry, the matrix elements of the local currents, both 
leading and subleading, must be written in terms of the effective states \M(v)). However, 
these states are not eigenstates of the Hamiltonian, once Ok and Oq are included in the 
Lagrangian. Hence we must allow for the possibility that if an effective state \M(v)) is 
created at time t = — oo, then Ok or Oq could act on the state before its decay at t = 0. 
This possibility is accounted for by including time-ordered products in which Ok or Oq is 
inserted along the incoming or outgoing heavy quark line. If we are keeping terms of order 
1/toq, only one insertion of Ok or Oq needs to be included. The time-ordered products are 
of the form JT3| 



(D^(v)\cTb \B(v)) = . . . + JL (M'(v')\i J dyT {h c v ,rh b v , O k + 0%) \M{v)) 

+ (-M W / dyT{K,Th b v , O k + O^} \M{v)) , (3.22) 
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where the ellipses refer to the current corrections computed earlier. The evaluation of 
the matrix elements of the time-ordered products will lead to still more nonperturbative 
functions like F(v,v') and G a (v,v'). 



C. Corrections at zero recoil 

It is straightforward, but not very illuminating, to expand all of the new nonperturba- 
tive functions which arise at order I/tuq in terms of scalar form factors. In the end, the 
corrections may be parameterized in terms of four functions of the velocity transfer w, and 
a single nonperturbative parameter A, all proportional to the mass scale Aqcd- The new 
parameter has a simple interpretation as the "energy" of the light degrees of freedom, and 
is given by 

A = lim (rriB — rrib) ■ (3.23) 

nib— >oo 

Instead of a general treatment, however, we will consider the I/ttlq corrections at the 
zero recoil point w — 1. This is clearly the most important case, because it is at this point 
that the nonperturbative matrix elements are absolutely normalized in the heavy quark 
limit. What happens to this normalization condition when 1/wiq corrections are included? 

Let us study the corrections to the current in detail. They are described by the nonper- 
turbative function G a (v,v'). At v — v', G a (v,v) may be expanded as 

G a (v, v) = G lVa + G 2la + G 3 vJ + G^J . (3.24) 

But G a (v,v) is subject to the same constraint as F(v,v'), 

G a (v,v) = P_G a (v,v)P- = (d - G 2 - G 3 + G A )v a P- = Gv a P_ , (3.25) 

and it, too, is equivalent to a Dirac scalar (the same is not true of the general function 
G a (v,v')). Now consider the matrix element where we take T a = v a . Then we have 



(M'(v)\h v/ ivDh v \M(v)) = - [M M M M ,} 1/2 Tr \m\v) v a M(v) 



Gv 



a 



- G [M M M M ,} 1/2 Tr [M'(v)M(v)] . (3.26) 
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But this matrix element vanishes by the classical equation of motion in the effective theory, 

v-Dh v (x)=0. (3.27) 

Hence G = = G a (v,v). There are no 1/toq corrections from the current to the normaliza- 
tion condition at zero recoil Ji"3"| , |2"U . 



The same is true of insertions of the corrections Ok and Oq to the Lagrangian: their 
contribution vanishes at w = 1. To show this requires the imposition of the conservation of 
the b number current at order l/m&, much as we derived the normalization of the Isgur-Wise 



function at leading order. The complete derivation is found in the literature [13]. 

In the end, we have an extremely exciting result, known as Luke's Theorem. There are 
no corrections at zero recoil to the hadronic matrix elements responsible for the semilep- 
tonic decays B — > D tv and B — > D* £v. The leading power corrections to the normalization 
of zero recoil matrix elements are only of order l/m? c . Given that Aqcd/^c ~ 30% and 
AQ CD /m^ ~ 10%, the implication is that the leading order predictions at w — 1 are consid- 
erably more accurate than one might have expected. In addition, away from zero recoil the 
l/m c corrections must be suppressed at least by (w — 1). 

On closer inspection, this result is more interesting for B — > D* Ev than for B — > D Iv. 
This is because the leading order matrix element for B — > D iv vanishes kinematically at 
zero recoil for a massless lepton in the final state. Hence, in this case the l/m c corrections 



are not suppressed as a fractional correction to the lowest order term [21 . 



D. Extraction of \V c b\ from B — » D* Iv 

An immediate application of these results is the extraction of \V c b\ from the exclusive 
decay B — > D* Iv. This process is mediated by the weak operator Ob c Ql-1] ), whose matrix 
element factorizes as 

(D* tv\ O bc \B) = (D*\ cy(l - 7 5 )& \B) (£v\ l llx {l - 7 > |0) . (3.28) 
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The leptonic matrix element may be computed perturbatively, while we treat the hadronic 
matrix element in the heavy quark expansion. The result is a differential decay rate of the 



form 22 



* 1 ' ( ' h ^V cb \ 2 (m B - m D *) 2 m 3 D *(w + i)Vw 2 - 1 



dw 487T 3 

x 



^ Aw m 2 B — 2wmbm£)* + m 2 D , 



F 2 {w) . (3.29) 



w + 1 (rag — m D *) 2 
All of the HQET analysis goes into the factor F(w), which has an expansion 

F(w) = £(w) + (radiative corrections) + (power corrections) . (3.30) 

We extract \V c b\ by studying the differential decay rate near w — 1, where the hadronic 
matrix elements are known. Of course, this requires extrapolation of the experimental data, 
since the rate vanishes kinematically at w = 1. For massless leptons, only the matrix element 
(D*\ cypj b \B) of the axial current contributes at this point. The analysis of this quantity 
in the HQET yields an expansion of the form 



^(1) = VA 





H 1 1" 0l/ m 2 + 

m c nib 



(3.31) 



The correction 5i/ m a, which contains terms proportional to l/m 2 , l/m 2 and l/m c mf,, is 
intrinsically nonperturbative. It has been estimated from a variety of models to be small 
and negative 0^3, 

8 1/m 2 w -0.055 ± 0.035 . (3.32) 

Note that the model dependence in the result has been relegated to the estimation of the sub- 
subleading terms. The radiative correction t]a has now been computed to two loops ]25 , 26 



tia = 0.960 i 0.007 . (3.33) 
The result is a value for F(l) with errors at the level of 5%, 

F(l) = 0.91 ±0.04. (3.34) 
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This is the theory error which the experimental determination of \V c b\ will inherit. It is 
dominated by the uncertainty in the nonperturbative corrections, and it is difficult to see 
how this can be improved much in the future. 

All that is left experimentally is to extrapolate the data to w = 1 and extract 

lim ; 1 ^- . (3.35) 
w ~*i \/w — 1 dw 

Once the kinematic factors in Eq. (|3.29|) have been included, this amounts to a direct 
measurement of the combination |V c fe|F(l). Four experiments recently have reported results 



for this quantity [|27j --|3C 



ALEPH : 


(31.4 ±2.3 ±2.5) x KT 3 




ARGUS : 


(38.8 ±4.3 ±2.5) x KT 3 




CLEO : 


(35.1 ±1.9 ±2.0) x 1(T 3 




DELPHI : 


(35.0 ± 1.9 ±2.3) x 10~ 3 . 


(3.36) 



Neubert |3l] has combined these measurements to find a "world average" of (35.3 ± 1.8) x 
10~ 3 , from which we extract 

\V cb \ = (38.8 ± 2.0 cxp ± 1.6th) x 10- 3 . (3.37) 

This value of \V c b\ has almost no dependence on hadronic models. In contrast to model-based 
"measurements", here the theoretical error is meaningful, in that it is based on a systematic 
expansion in small quantities. 



IV. INCLUSIVE B DECAYS 

An exclusive semileptonic B decay, such as B — ► D £v, is one in which the final hadronic 
state is fully reconstructed. An inclusive decay, by contrast, is one in which only certain 
kinematic features, and perhaps the flavor, of the hadron are known. In this case, we need 
a theoretical analysis in which we sum over all possible hadronic final states allowed by the 
kinematics. Fortunately, this is possible within the structure of the HQET. 
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As in the case of exclusive decays, the key theme is the separation of short distance 
physics, associated with the heavy quark, from long distance physics, associated with the 
light degrees of freedom. We will also rely on heavy quark spin and flavor symmetry. 
However, the new ingredient will be the idea of "parton-hadron duality", which, as we 
will see, also relies on the heavy quark limit m b 3> Aqcd- 



A. The inclusive decay B — > X c iu 

Let us consider the inclusive decay 

B(p B )^X c (p x )e(p e )is(p v ), (4.1) 

where all that is known about the state X c is its energy and momentum, and the fact that 
it contains a charm quark. This decay is mediated by the weak operator Ob c - It is easy to 
generalize our discussion to inclusive decays of other heavy quarks, such as b — > uiv and 
c — ► slu, by replacing (9 fcc with the appropriate weak operator. 

The treatment of exclusive decays required both the b and c quarks to be heavy. For 
inclusive decays we can relax this condition on the c quark, requiring only m& ^> Aqcd- 
What does the weak decay of the b, at time t = 0, look like to the light degrees of freedom? 
For t < 0, there is a heavy hadron composed of a point-like color source and light quarks and 
gluons. At t — 0, the point source disappears, releasing both its color and a large amount of 
energy into the hadronic environment. Eventually, for t > 0, this new collection of strongly 
interacting particles will materialize as a set of physical hadrons. The probability of this 
hadronization is unity; there is no interference between the hadronization process and the 
heavy quark decay. There are subleading effects in powers of AQ CD /m b , but they do not 
alter the probability of hadronization. Rather, they reflect the fact that the b quark is not 
exactly a static source of color: it has a small nonrelativistic kinetic energy and it carries a 
spin, both of which affect the kinematic properties of its decay. 

As in the case of exclusive decays, we will compute the inclusive semileptonic width 
F(B — > X c tv) as a double expansion in a s (rrib) and AQ CD /m fe [|ll^ ,j32}-|34l| . The expansion 
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in a s (mb) reflects the applicability of perturbative QCD to the short distance part of the 
process. The heavy quark expansion will be continued to relative order 1/ml, as there is an 
analogue of Luke's Theorem which eliminates power corrections to the rate of order l/m*,. 
The 1/ml corrections will be written in terms of three nonperturbative parameters. The 
first, A, is defined in Eq. ( |3.23| ). It is essentially the mass of the light degrees of freedom in 



the heavy hadron, but we will see that it is plagued by an ambiguity of order Aqcd in the 
definition of the b quark mass. The other two parameters are the expectation values in the 



B meson of the leading corrections Ok and Oq to £hqet- They are defined as |L7 



Ai = ^—{B\ O k \B) 
2m B 

\2 = -^-(B\O g \B), (4.2) 
bras 

where we take the usual relativistic normalization of the states. Hence, Ai may be thought 
of roughly as the negative of the h quark kinetic energy, and A2 as the energy of its hyperfine 
interaction with the light degrees of freedom. 

Now let us outline the computation. The inclusive decay involves a sum over all possible 
final states, which is actually a sum over exclusive modes (such as D, D*, Dtt, . . .), followed 
by a phase space integral for each mode. We write 

T{B^XJv)=Y J I d[P-S.] \(X c £u\O bc \B)\ 2 . (4.3) 

There is an Optical Theorem for QCD, which follows from the analyticity of the scattering 
matrix as a function of the momenta of the asymptotic states. Its content is that a transition 
rate is proportional to the imaginary part of the forward scattering amplitude with two 
insertions of the transition operator, 

Y(B^X c iu) = -2Imz J dx e ikx (B\T {ol{x),O bc {0)} \B) = 2lmT. (4.4) 

In what follows, we will write the time-ordered product T{0\ c , Ob c } as a series of local 
operators, using the Operator Product Expansion. As we will see, the applicability of this 
expansion, and its computation in perturbation theory, will rest on the limit m b ^> Aqcd- 
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We will then use this limit again to expand the matrix elements of these local operators in 
the HQET. 

The first step is to factorize the integration over the lepton momenta, which can be 
performed explicitly. Written as a product of currents, O bc takes the form 

O bc =^J b U £lx , (4.5) 

where 

Jl = &f(l - 7 5 )& 

J$ = tf{l—f)v. (4.6) 

Then T can be decomposed as an integral over the total momentum = transferred 
to the leptons, 

T = \G 2 F \V cb \ 2 J dqT^(q)L^(q) . (4.7) 

Here the lepton tensor is 



L, v (q)= /d[P.S.] (0\Jl\£u) (iu\J eu \0) 

= ^ (q»Qu - q 2 g»v) , (4.8) 

and the hadron tensor is 

T^(q) = -ij dxe^(B\T{j£(x),Ji; c (0)} \B) . (4.9) 

We will need the imaginary part, Im T^ v '. Where is it nonvanishing? In quantum field theory, 
a scattering amplitude develops an imaginary part when there can be a real intermediate 
state, that is, the intermediate particles can all go on their mass shell. Whether this is 
possible, of course, depends on the kinematics of the external states. 



In this case, there are two avenues for creating a physical intermediate state The 
first is to act on the external state \B) with the transition current J bc . The state which is 
created has no net b number and a single charm quark; the simplest possibility is the decay 
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process b — *> c. The momentum of the intermediate state is px = Pb — q~, the condition that 
it could be on mass shell is simply 

Px = (PB-q) 2 >m 2 D - (4-10) 

If we define scaled variables 

p B = m B v fl , qf M = qf J -/m B , m D = m D /m B , ( 4 -H) 
this condition becomes 

vq<^(l + q 2 -m 2 D ) . (4.12) 

Another possibility is to act on \B) with the conjugate operator J{£ ■ This operation would 
produce an intermediate state with two b quarks and one c. For this state to be on shell, 
the momentum transfer has to satisfy 

Px = (Pb + qf > {2m B + m D f , (4.13) 

that is, 

v-q> i (3 -<f + 4m D + m^) . (4.14) 

The physical intermediate states are shown as cuts in the v ■ q plane in Fig. ^. Also shown 
is the contour corresponding to the phase space integration over the lepton momentum q. 
For physical (massless) leptons which are the product of a heavy quark decay, this integral 
runs over the top of the lower cut, for the range 

Iq^ + ie < v -q < ^ (l + <f - m 2 D ) + ie . (4.15) 



2 

As indicated by the dotted line, we can continue this contour around the end of the cut and 
back along the bottom, to v ■ q — y/q 1 — ie. Since T^ u {y ■ q *) = —T^ v (v ■ q) for real q 2 , we 
compensate for extending the contour by dividing the new integral by two. 

We now encounter our central problem. The integral over v ■ q runs over physical interme- 
diate hadron states, which are color neutral bound states of quarks and gluons. Hence the 
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v ■ q 


integration contour 

I (extended) 








physical cut / 




unphysical cuts 



FIG. 4. Analytic structure of T^ u in the complex v ■ q plane, for fixed, real q 2 . The integration 
contour is over the "physical cut", corresponding to real decay into leptons. The unphysical cuts 
correspond to other processes. 

integrand depends intimately on the details of QCD at long distances, which is intrinsically 
nonperturbative. A perturbative calculation of T^, which is all we have at our disposal, 
would appear to be of no use. 

The solution is to deform the contour away from the cut, into the complex v ■ q plane, 
as shown in Fig. |5|. Since the scale of momenta is set by m&, the contour is now a distance 



~ rrib away from the resonances |32| . Since m& ^> Aqcd, it is reasonable to hope that 
a perturbative treatment in this region is valid. Essentially, we are saved because we do 
not need to know T^ u {q) for every value of q, just suitable integrals of T^ v . That we can 
use such arguments to compute perturbatively the average value of a hadronic quantity, 
where at each point the quantity depends on nonperturbative physics, is known as (global) 
parton-hadron duality. 

Parton-hadron duality has the status of being somewhat more than an assumption, since 
it is known to hold in QCD in the limit mj, — > oo, but somewhat less than an approximation, 
since it is not known how to compute systematically the leading corrections to it. In any 
case, the limit m& ^> Aqcd plays a crucial role here. By deforming the integration contour a 
distance ~ away from the resonance regime, we find the correspondence in QCD of our 
earlier intuitive stateprobabilityment: the probability of the decay products materializing as 
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FIG. 5. The deformation of the integration contour into the complex v ■ q plane. 



Q q 

Pb-q 





FIG. 6. The operator product expansion at tree level. 

physical hadrons is unity, independent of the kinematics of the short distance process. The 
local redistribution of probability in phase space due to the presence of hadronic resonances 
is irrelevant to the total decay. Finally, we should note that since we do not have control 
over the corrections to local duality, it might work better in some processes than in others, 
for reasons that need not be apparent from within the calculation. Hence one must be par- 
ticularly wary of drawing dramatic conclusions from any surprising results of these inclusive 
calculations |35| . 

Let us perform the operator product expansion at tree level, and for decay kinematics. 
The Feynman diagram is given in Fig. |6], which yields the expression 

i> h - 4 + m c 



{Pb — q) 2 —ml + ie 



(4.16) 



We now write 



= q^/m b 
rh c = m c /m b 
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b(x) 



v - x h v {x)+0(l/m b ) 



(4.17) 



and expand in powers of l/m b . Since the operator product expansion is in terms of the 
effective field h v , a factor of corresponds to an insertion of the covariant derivative iD^. 
Operator ordering ambiguities are to be resolved by considering graphs with external gluon 
fields. 

As an example of the procedure, let us expand the propagator to order \jm\. (There 
are also corrections to the currents at this order, which are included in a full calculation.) 
It is convenient to define the scaled hadronic invariant mass, 



{Pb - qf '/ml = 1 - 2v ■ q + q 2 . 



(4.18) 



Then we find a contribution to T^ v of the form 



TV" = — h vl ^f - D7*(l - 7 5 ) 



1 



2k ■ q — q 2 



h v + .... (4.19) 



s — m 2 + ie (s — m 2 + ie) 2 
From this expression we can read off the operators which appear in the operator product 
expansion. Since 

1 



Im 



6^(3 - m 2 c ) , 



(4.20) 



(s — m 2 c + ie) 7 

we see that the effect of taking the imaginary part in each term is to put the charm quark 
on its mass shell. The leading term is a quark bilinear, 

1 



m b 



h v Y^-ih u (i-r)K 



(4.21) 



It is straightforward to compute its matrix element in the HQET using the trace formalism, 
(B\ KY(1> - |)t' / (1 - 7 5 K, \B) = 2m B {2v^v u - gT - ^<f - v v q? + v ■ q + . . .) . 

(4.22) 

The ellipses denote terms of order \jm\\ there are no corrections of order l/m b , by Luke's 
Theorem. Finally, we contract the tensor with and perform the phase space integra- 
tion ( |4.7|) . In the end, the result is the same as we would have gotten directly by computing 
free quark decay. 
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Of course, if we only intended to reproduce the free quark decay result, we would never 
have introduced so much new formalism. The value of the HQET framework is that it 
allows us to go beyond leading order and compute the next terms in the series in l/m%. For 



example, consider the operators induced by the expansion of the propagator ( |4.19| ). The 
correction of order 1 / m^ comes from the operator 



1 



mi (s — mi + ie) 



K-fy-fi-fil-^q-iDhv. 



However, the matrix element of this operator is of the form 

(B\h v T a (v,q)iD a h v \B), 



(4.23) 



(4.24) 



which, as we have seen, vanishes by the classical equation of motion. In fact, since all l/m^ 
corrections, from any source, have a single covariant derivative, they all vanish in the same 
way. This is the analogue of Luke's Theorem for inclusive decays |52]. The correction of 
order \jm\ in Eq. ( [4.1 9| ) is 

1 1 



mj, + ie) 



(4.25) 



mi [s 

The matrix element of this operator is related by the heavy quark symmetry to Ai, the 
expectation value of Ok- The full expansion of T^ u also induces operators with explicit 
factors of the gluon field, whose matrix elements are related to A2. 

We now present the result for the inclusive semileptonic decay rate, up to order in 
the heavy quark expansion, and with the complete radiative correction of order a s . We also 
include that part of the two loop correction which is proportional to (3 Q a 2 s . Since (3 Q ~ 9, 
perhaps this term dominates the two loop result. In any case, it is interesting for other 
reasons, as we will see below. 

Let us first consider the decay B — > X u 
m u = 0. We find 



for which the decay rate simplifies since 



G 2 F \V u b\ 2 

1927T 3 



ml 





/25 


2vr 2 \ 


1 + 


l"6 


3 ) 



71 



(2.98/?o + C u ) 



'a s (m b y 



71 



Ai - 9A 2 
2m\ 



+ ... 



+ ... 

(4.26) 



35 



When we include the charm mass, it is convenient to write the unknown quark masses in 
terms of the measured meson masses and the parameters of the HQET. In terms of the spin 
averaged mass Tub = {itlb + 3ms* )/4, we have 



777,6 = 777b — A + 



_Ai_ 
2m f 



+ 



(4.27) 



and analogously for m c . We then find |3l|j34|j38|j39 
T(B 



n%\\r |2 

UF]Vcbl m%x 0.369 



192tt3 

-1.65 A fi_ .87^~ 



1.54^- (1.43/3o + C C ) f^l 



7T 



0.95- 



m 



B 



3.18^ + 0.02- 



TT 

A2 

2 

ml 



+ . . . 



(4.28) 



All the coefficients which appear in this expression are known functions of mu/rfiB, and 
are evaluated at the physical point = 0.372. In both B — > X u tv and B —* X c £u, 

the power corrections proportional to Ai and A2 are numerically small, at the level of a few 
percent. 



B. Renormalons and the pole mass 

The inclusive decay rate depends on the heavy quark mass m&, either explicitly, as in 
Eq. (|4.26 ), or implicitly through A, as in Eq. (|4.28| ). At tree level, ra& is just the coefficient 



of the bb term in the QCD Lagrangian, but beyond that we are faced with the question of 
what exactly we mean by m&. Should we take an MS mass, such as m^Wi,)? Or should we 
take the pole mass ml° le , or maybe some other quantity? The various prescriptions for m& 
can vary by hundreds of MeV, and, since the total rate is proportional to ml, the question 
is of practical importance if we hope to make accurate phenomenological predictions. 

At a fixed order in QCD perturbation theory, the answer is clear. The heavy quark 
masses which appear come from poles in quark propagators, so we should take mf ole (and 
mP ole ). This is also the prescription for the mass which cancels out the on-shell part of the 
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heavy quark field in the construction of £hqet- Hence the difference of heavy quark pole 
masses is known quite well, 

m r-<-=(M B -A + ^ + ...)-(^-A + A^ + ...) 

= 3.34GeV + O(A?, 0D /m|). (4.29) 

Since T(B — > X c tv) depends approximately as ml(m b — m c ) 3 , the uncertainty due to quark 
mass dependence is reduced. 

The problem, of course, is that there is no sensible nonperturbative definition of m b ole , 
since due to confinement there is no actual pole in the quark propagator. Hence a direct 
experimental determination of a value for mf ole to insert into the theoretical expressions 
( |4.26| ) and ( }4.28|) is not possible. How, then, can we do phenomenology? 

One approach would be to define m b ole to be the pole mass as computed in perturbation 
theory, truncate at some order, and then estimate the theoretical error from the uncomputed 
higher order terms. However, it turns out that even within perturbation theory the concept 
of a quark pole mass is ambiguous. Consider a particular class of diagrams which contribute 
to m% ole , shown in Fig. 0. The perturbation theory is developed as an expansion in the small 
parameter a s (rrib), so we hope that it will be well behaved. Each of the bubbles represents 
an insertion of the gluon self-energy, which is proportional at lowest order to a s (m b )(3 . Of 
course, the infinite sum of the graphs in Fig. [7| can be absorbed into the one loop graph, with 
a compensating change in the coupling from a s (m b ) to a s (q), where q is the loop momentum. 
The result is an expansion for m^ ole of the form 

1 + aia s (m b ) + (a 2 /3 + b 2 )a 2 s (m b ) + (a 3 (3 2 + b 3 (3 + c 3 )a^(m 6 ) + ... . 

(4.30) 

The graphs in Fig. [?| contribute the terms proportional to a™ + (m b ) (3$ . Since j3 ~ 9 these 
terms are "intrinsically" larger than ones with fewer powers of (3q, and we might hope that 
their sum approximates the full series. However, it is important to realize that the only limit 
of QCD in which such terms actually dominate is that of large number of light quark flavors, 
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m P° e = m b (m b ) 



+ A + ... 

FIG. 7. The radiative corrections to m-b° le °^ or der a s {mi ) ) n+l (5q . 

in which case the sign of (3q is opposite to that of QCD. Although this is a physical limit of 
an abelian theory, we are certainly not close to that limit here. The ansatz of keeping only 
the terms proportional to a™ +1 (rrib)f3o is known as "naive nonabelianization" (NNA) [fflj|. 

What is most interesting about the series of terms shown in Fig. [7|, a n aC ( m 6)/3o _1 > * s 
that it does not converge. Already in the graphs kept in the NNA ansatz, we are sensitive 
to the fact that QCD is an asymptotic, rather than a convergent expansion. For large n the 
coefficients a n diverge as n\, much stronger than any convergence due to the powers a™ (rat). 
The series can only be made meaningful if this divergence is subtracted. As with many 
subtraction prescriptions, there is a residual finite ambiguity^ This ambiguity, known as an 



"infrared renormalon", leads to an ambiguity in the pole mass of order [41-43 



5m p b ole ~ 100 MeV . (4.31) 

By the definition ( |3.23| ), A also inherits this ambiguity. 

The expressions ( J4.26|) and ( |4.28| ) are plagued by two problems. The first is the renor- 



malon ambiguity in m^ ole and A. The second is that the perturbative expansion for the rate 
r is itself divergent, and also has an infrared renormalon. In the expansion 

r = r a' n a^{m^)^Q ~ 1 + (power corrections) , (4.32) 



*In a formal treatment, this ambiguity arises from a choice of contour in the Borel plane. 
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the coefficients a' n also diverge as n\. However, it turns out that these two problems actually 
cure each other, because the infrared renormalons in ml° le and in the perturbation series 



for T cancel |42|{l4j • We can exploit this cancelation to improve the predictive power of the 
theoretical computation of the rate. Without this improvement, the infrared renormalons 
render the expressions ( |4.26| ) and Q4.28| ) of dubious phenomenological utility. 

The most reliable approach, theoretically, is to eliminate m£ ole or A explicitly from the 
rate by computing and measuring another quantity which also depends on it. For example, 
let us consider the charmless decay rate T(B — > X u iv) and the average invariant mass 
(sh) of the hadrons produced in the decay. Each of these expressions suffers from a poorly 
behaved perturbation series in the NNA approximation. Ignoring terms of relative order 
I /ml and writing the rate in terms of A instead of m£ ole , we find to five loop order [41], 



G^|Kb| 2 5 

T = — rrir, 

192vr 3 3 



1 _ 2.41^ - 2.98 ( ^ ) p ~ 4-43 ( ^ 



01 



7.67 



15.7 ('^V $ + ...- 5— + 
7r / m B 



G 2 F \V U b\ 

1927T 3 



■ m 



B 



1 - 0.061 - 0.120 - 0.107 - 0.111 - 0.136 + . . . 
- 5A/m B + ... 



(4.33) 



for a s {m b ) = 0.21 and /3 = 9. As we see, not only does the perturbation series fail to 
converge, it does not even have an apparent smallest term, where one should truncate to 
minimize the error of the asymptotic series. The series for (sh) exhibits a similar behav- 
ior gp], 

2 /™ \ 3 /™ \ 4 



(s H ) = m 



B 



0.20^ + 0.35 i 3 + 0.64 (^') /3 2 + 1.29 (^) 0, 

+ 2.96 + + 1-5- 

V 7r / 10 m B 



+ ... 



m 



B 



0.0135 + 0.0141 + 0.0156 + 0.0189 + 0.0261 + . . . - 7A/10m B + . . . . (4.34) 



However, the situation improves dramatically if we eliminate A and write T directly in terms 
of (s H ), 
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G 2 F \V U b 



2 



1 = — m 



,5 
B 



1 - 7.14^- - 0.064 - 0.020 - 0.0002 - 0.022 - 0.047 + . . . 



(4.35) 



1927T 3 

By truncating this series at its smallest term, 0.0002, we obtain a new expression in which 
the theoretical errors are under control. The price is that we must now measure a second 
quantity, (sh), in the same decay. 

An analogous approach works for inclusive decays to charm, in which case the new 
quantity which we must measure is (sh — Let us apply this analysis to the extraction 
of V c b from the inclusive decay B — > X c tv. Eliminating A as before, we find El 



Y[B -> XJtv) = ^^r™l x 0.369 



a s {m b ) a 2 (m b ) 



1 - 1.17 v ' - 0.74/3, 



AT o 
7T 7T Z 



,(s H -m 2 D ) 



-7.17 ' 9 + 



(4.36) 

plus small terms of order Xi^/rri^. The current data on (s# — mb) are still inconclusive, 
but they will improve in the future. As an illustration of the utility of the method, however, 
let us consider the physical limit (sjj — m 2 D ) > 0. In that case, we already can obtain the 
bound 

/ \ -1/2 

V cb > 0.037 x — , (4.37) 
y 1.54ps / 

where we have used the observed B semileptonic branching ratio of 10.7%. As the data on 
(sh — ttL 2 d) improve, so will this bound. Note that it is entirely consistent with the value 
( [3.371 ) obtained from exclusive decays. 

An alternative approach is to express the width V in terms of the running mass mb(mb) 



instead of another inclusive observable | 42| , |45 l. Since the MS mass is a short distance quan- 



tity, this also eliminates the infrared renormalon, which is associated with long distance 
physics. However, from a phenomenological point of view, it begs the question of how the 
running mass is to be determined from experiment. Possibilities include quarkonium spec- 
troscopy, QCD sum rules, and lattice calculations, but in all of these cases it is unclear both 
how to determine reliably the accuracy of the method, and how to deal with renormalon 
ambiguities in a manner that is consistent with their treatment in the calculation of T. With 
additional progress, however, such an approach might eventually prove fruitful. 

40 



V. HEAVY HADRON PRODUCTION VIA FRAGMENTATION 



So far we have applied the idea of heavy quark symmetry to heavy hadron spectroscopy 
and to heavy hadron decay. We can also apply it to heavy hadron production via frag- 
mentation. Although there is not yet available the wealth of detailed data for production 
processes that there is for spectroscopy and decays, the study of heavy quark fragmentation 
is an elegant application of the heavy quark limit in a new regime ||46|| . 

A. The physics of heavy quark fragmentation 

The production of a heavy hadron proceeds in two steps. First, the heavy quark itself 
must be created; because of its large mass, this process takes place over a time scale tq 
which is very short, tq < 1/itlq. Second, some light degrees of freedom assemble themselves 
about the heavy quark to make a color neutral heavy hadron, a process which involves 
nonperturbative strong interactions and typically takes much longer, t f ~ 1/Aq CD ^> tq. If 
the heavy quark is produced with a large velocity in the center of mass frame, and if there is 
plenty of available energy, then production of these light degrees of freedom will be local in 
phase space and independent of the initial state. This is the fragmentation regime. We will 
see that heavy quark symmetry simplifies the description of heavy hadron production via 
fragmentation, because, as before, it allows us to separate certain properties of the heavy 
quark from those of the light degrees of freedom. This is particularly important in the 
production of excited heavy hadrons, for which the behavior of the spin of the light degrees 
of freedom can be quite interesting. 

Heavy quark symmetry simplifies the picture of fragmentation in several ways. First, 
since the fragmentation process involves the exchange only of soft momenta, the velocity 

of the hadron Hq which is produced is the same as the velocity of the heavy quark 
Q. Second, the angular momenta Sq and Jt of the heavy quark and the light degrees of 
freedom are well defined and can be followed individually through the processes of production 
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and decay. Third, the typical precession time ts of Sq and Jt about each other is of 
order r s ~ ""^qMqcd 1/Aq CD , since the precession is induced by the chromomagnetic 
interaction. Hence rs t>, and during the fragmentation process Sq can essentially be 
taken to be frozen in direction. 

Suppose that Q fragments initially to a state with light degrees of freedom with spin 

Jg. This could correspond to two possible hadronic states: Hq, with spin J = Jg — |, or 

1 (*) 

Hq, with spin J = Jc + ±. Let the states Hq have lifetime r. How should we treat the 
interaction of Sq with Jp. 

The answer depends critically on the relationship of r to the precession time ts- While 
the order of magnitude of r s is essentially fixed at ?71q/Aq CD , the lifetime r can vary widely, 

depending on the channel in which Hq^ primarily decays. Of course, if Hq^ decays weakly, 

(*) 

then t 3> ts- But if Hq decays strongly, then r depends critically on the phase space 
available for the transition. In the absence of phase space factors, r ~ 1/Aq C d <C Ts- But 
since the pion, the lightest hadron, is not massless, it is possible to have a situation where 
an allowed strong decay is either forbidden by phase space (as with B* — > B + n), or severely 
suppressed (as with D* — > D + it). 

The treatment of a given hadron doublet depends on the relative size of r and ts- Let us 
consider the two extreme possibilities. The first corresponds to a strong decay with plenty 
of phase space, so r s > r. Here Hq and Hq are formed and then decay before the angular 
momenta Sq and Jg have a chance to interact. The fact that there is no precession means 
that there is no change in the orientation of Sq and Jg. If, for example, either spin were 

polarized before or during the fragmentation process, this polarization would be undisturbed 

(*) 

by passing through the resonance Hq' . The strong decay of the heavy hadron then would 
reflect the orientation of the light degrees as they were produced. 

Note that it is the very same spin exchange interaction which is inhibited here which is 
responsible for the splitting A H between Hq and H Q . Hence, under these conditions the 
resonances are almost completely overlapping, with widths T = 1/t satisfying T ^> A#. This 
is another consequence of the effective decoupling of Sq and Je, which are independent good 

42 



quantum numbers of the resonances. It also provides a convenient criterion for determining 
directly whether we are in the regime r s ^> r: we cannot distinguish separate resonances 
for H Q and H Q . 

The second possibility is the opposite extreme, t s <C r. This corresponds to heavy 
hadrons which decay weakly or electromagnetically, or to strong decays which are severely 
suppressed by phase space. Here the spins Sq and Jp have plenty of time to interact, 
precessing about each other many times before Hq and Hq decay. There is at least a partial 
degradation of any initial polarization of Sq, as well as a degradation of any information 
about the fragmentation process which may be carried by the light degrees of freedom. The 
signature of this situation is that the states Hq and Hq are well separated resonances, since 
the chromomagnetic interactions have ample opportunity to produce a hyperfine splitting 
much larger than the width, Ah ^> T. In contrast with the first case, here the heavy and 
light spins are resolved into states of definite total spin J = \ Jp ± \ \ . This resolution destroys 
information about the individual angular momenta Sq and Jp, and the decays of Hq and 
Hq only partially reflect the features of the production process. 

B. Production and decay of Di and D\ 

We will illustrate these ideas with two examples. The first concerns the production 
and decay of the excited charmed mesons D\ and D\. These are states with light degrees 
of freedom with angular momentum and parity J[ = | + ; in the quark model, this would 
correspond to a light antiquark with a unit of orbital angular momentum. The physical states 
are the D u with mass 2423 MeV and width 18MeV, and the D* 2 , with mass 2458 MeV and 
width 21 MeV. The splitting between them is 35 MeV; although this is not much larger than 
the individual widths, for simplicity we will treat them in this limiting case r ^> Ts- 

Let us review the sequence of events. First, the charm quark is created in some hard 
interaction, with a time scale r c ~ l/m c . Second, light degrees of freedom with Jp = | 
are created in a fragmentation process, over a time tf ~ 1/Aqcd, forming a color neutral 
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charmed hadron. Third, the angular momenta S c and Je precess about each other, resolving 
the individual D\ and D\ states, which requires a time t$ ~ m c /AqcD- Then, finally, after 
a time r greater than m c /AQCD, the D% or D\ emits a pion and decays to a D or D*. One 
question we might ask about this process is, what is the alignment of the light degrees of 
freedom produced in the fragmentation? Of course, the answer depends on nonperturbative 
QCD and is difficult to calculate from first principles. Still, we can explore the question 
experimentally. 

Since the total angular momentum of the light degrees of freedom is Jt = |, they can 
be produced with helicity h — ±| or h — ±| along the fragmentation axis. While parity 
invariance of the strong interactions requires that the probabilities for helicities h and —h 
are identical, the relative production of light degrees of freedom with \h\ = | versus \h\ = \ 
is determined in some complicated way by nonperturbative dynamics. Let the quantity W3/2 



denote the probability that \h\- ' 



2 ■ 

3\ 1 TD(U _ 3n 



w 3/2 = P(h=$)+P{h = -l). (5.1) 

Then 1 — u> 3 / 2 is the probability that \h\ = \. Completely isotropic production corresponds 
to W3/2 = \- The observable W3/2 is a new nonperturbative parameter of QCD, which is well 
defined only in the heavy quark limit. 

This new parameter can be measured directly in the strong decay of the D\ or D\. For 
example, consider the angular distribution of the pion with respect to the fragmentation 
axis in the decay D\ — > D + 71. This is a decay of the light degrees of freedom in the excited 
hadron, so it will depend on their initial orientation (that is, on W3/2) and on the details of 
the precession of Jt around Sq during the lifetime of the D\. Following the direction of 



through the sequence of fragmentation, precession and decay, we find the distribution [|| 

1 dr x 



Tdcosfl 4 



1 + 3 cos 2 9 - 6w 3/2 (cos 2 9 - |)J . (5.2) 

This distribution is isotropic only when W3/2 = |, that is, when the light degrees of freedom 
are produced isotropically in the fragmentation process. Similar distributions are found in 
the decays D\ — > D* + n and D 1 — > D* + 7r. 
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A fit of ARGUS data |47j to the expression (|5.2|) seems to indicate that a small value 
of w 3 /2 is preferred; while the errors are large, we find that w 3 / 2 < 0.24 at the 90% confi- 
dence level. It would be nice to confirm this result with a sharper measurement, and not 
only for the charmed mesons but in the bottom system as well. Since W3/2 is intrinsically 
nonperturbative, we do not have any real theoretical understanding of why it should be 
small, although one can construct a qualitative understanding based on a string picture of 



fragmentation [4S 



C. Polarization of A b at SLC/LEP 

The second example springs from a practical question: What is the polarization of A& 
baryons produced at the Z pole? This question is motivated by the fact that b quarks 
produced in the decay of the Z are 94% polarized left-handed. Since the A& is composed of 
a b quark and light degrees of freedom with Ji = 0, the orientation of a A& is identical to 
the orientation of the b quark inside it. Similarly, the b quark spin does not precess inside 
a Ab. Hence if a b quark produced at the Z fragments to a A&, then those baryons should 
inherit the left-handed polarization of the quarks and reveal it in their weak decay. 

Unfortunately, life is not that simple. Two recent measurements of A& polarization P(A&) 
from LEP are [H,[l 

DELPHI : 0.08 ^;g(stat.) +^(syst.) 
ALEPH : 0.26 ±g : |g(stat.) ±^(syst.) , 

both a long way from P(A&) = 0.94. The reason is that not all b quarks which wind up as 
A;, baryons get there directly. There is a competing process, in which they fragment first to 
the excited baryons and E£, which then decay to A& via pion emission. These new states 
have light degrees of freedom with J a — 1. If they have a lifetime r > rg, then the b quark 
will have time to precess about Ji and its polarization will be degraded. The result will be 
a net sample of A^'s with a polarization less than 94%, as is in fact observed. 
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In addition to the requirement that t > t$ for the any depolarization of A^'s by 
this mechanism depends on two unknown quantities. First, there is the production rate / of 
relative to A&. Isospin and spin counting enhance / by a factor of nine, while the mass 
splitting between £^ and A& suppresses it; studies based on the Lund Monte Carlo | 5"I[ 
indicate / pa 0.5 with a very large uncertainty. Of course, it would be preferable to measure 
/ directly. Second, there is the orientation of the spin Ji with respect to the fragmentation 
axis. This orientation, which is nonperturbative in origin, allows the helicities ft, = 1,0, — 1. 
By analogy with the treatment of the Di and let us define 

Wl = P(h=l)+P(h=-i). (5.3) 

In this case, isotropic production corresponds to w% = |. We may measure W\ from the 
angle of the pion with respect to the fragmentation axis in the decay — > A& + it |4(| , 

^t^-t = Hi + 3 cos 2 6 - | Wl (cos 2 9-1)}. (5.4) 
rdcosfc* 4 L 2 6 J 

It turns out that the decay £& — > A& + 7r is isotropic in cos# for any value of wi. 

The polarization retention of the Aj, may be computed in terms of / and w\. As before, 
we will consider the simpler situation in which the £& and the ££ do not overlap, sor > ts- 
Then the polarization of the observed A^'s is P(A&) = R(f,Wx)P(b), where P(6) = 94% is 
the initial polarization of the b quarks, and El 



l + l(l + 4 Wl )/ 

p(/, wi) = — iTJ ' ^ ^ 

Note that for / = (no s are produced), P(0, wi) = 1 and there is no depolarization. 
For the Lund value / = 0.5, R ranges between 0.70 and 0.85 for < W\ < 1. 

Can the very low measured values of P(A&) be accommodated by the present data on 
the ? The situation is still unclear, because not much is known about the E& and 
states. DELPHI ]52| has reported the observation of these states, but with masses which 
deviate significantly from those that one would expect based on heavy quark symmetry and 



the observed S c and S* [M - Along with the masses, DELPHI also reports w\ ~ and 
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1 < / < 2. If this is confirmed, then a polarization in the range -P(Afc) »s 40% — 50% is 
easy to accommodate. On the other hand, the CLEO analysis |kJ of the decay of the S* 
yields W\ = 0.71 ±0.13, consistent with isotropic fragmentation. Recall that by heavy quark 
symmetry, w\ measured in the charm and bottom systems must be the same, so this result 
is inconsistent with the report from DELPHI. Clearly, further measurements are needed to 
resolve this situation. 

VI. CONCLUDING REMARKS 

Unfortunately, we have had time in these lectures only to introduce a very few of the 
many applications of heavy quark symmetry and the HQET to the physics of heavy hadrons. 
Since its development less than ten years ago, it has become one of the basic tools of QCD 
phenomenology. Much of the popularity and utility of the HQET certainly come from its 
essential simplicity. The elementary observation that the physics of heavy hadrons can be 
divided into interactions characterized by short and long distances gives us immediately a 
clear and compelling intuition for the properties of heavy-light systems. The straightforward 
manipulations which lead to the HQET then allow this intuition to form the basis for a new 
systematic expansion of QCD. The deeper understanding of heavy hadrons which we thereby 
obtain will become increasingly important as the end of the second millennium approaches 
and the B Factory Era begins. 
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